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Image segmentation is often faced by low contrast, bad boundaries, and 
inhomogeneity that made it difficult to separate normal and abnormal tissue. 
Therefore, it takes long periodto read and diagnose brain tumor patients. The 
aim of this study was to applied hybrid methods to optimize segmentation 
process of magnetic resonance image of brain. In this study, we divide the 
brain tumor images with double density dual-tree complex wavelet 
transform (DDDTCWT), continued by convolutional neural network (CNN), 
and optimized by genetic algorithm (GA) with 48 combinations yielding 
excellent results. The F-1 score was 99.42%, with 913 images test data. The 
training images consist of 1397 normal magnetic resonance imaging (MRI) 
images and 302 tumor MRI images resized by 32 x32 pixels. The 
DDDTCWT transforms the input images into more detail than ordinary 
wavelet transforms, and the CNNs will recognize the pattern of the output 
images. Additionally, we applied the GA to optimize the weights and biases 
from the first layer of the CNNs layers. The parameters used for evaluating 
were dice similarity coefficient (DSC), positive present value (PPV), 
sensitivity, and accuracy. The result showed that the combination of 
DDDTCWT, CNN, and GA could be used to brain MRI images and it 
generated parameters value more that 95%. 
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1. INTRODUCTION 


Magnetic resonance imaging (MRI) is a radiological scanning technique that uses magnets, radio 
waves, and computers to produce images of body structures [1]. Various methods can be used for imaging a 
medical object, such as angiogram, brain scan, computerized tomography (CT)-scan, diffusion tensor 
imaging, functional magnetic resonance imaging (f(MRI) [2], MRI, magnetic resonance spectroscopy (MRS), 
positron emission tomography (PET), and Biopsy [3]. Some literature states that MRI is the best examination 
tool for its relatively safe radiation hazards [4] and the high accuracy rate. In an MRI examination, the patient 
is placed on a bed and inserted into a magnetic tube. A strong magnetic field will be formed and align the 
protons of the hydrogen atoms, which are then exposed to radio waves. The receiver would detect the signal 
on the MRI machine, then the computer would process the received information and produces an image. The 
images and resolution of MRI are presented in detail and can detect small changes in the body structures. In 
some procedures, a contrast material, such as gadolinium, is used to improve image accuracy. In addition, 
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MRI can demonstrate the presence of abnormal tissue with high resolution and good contrast [5]. 
Furthermore, MRI results have a higher sensitivity for detecting the presence or change in tumor size. A 
tumor is a disease that arises due to abnormal changes in body tissue cells that turn into tumor cells. When a 
brain tumor occurs, the excessive growth of unnecessary cells causes an increased volume and damage to 
other cells and interferes with the function of that part of the brain. Brain tumors can cause complications, 
including cerebral edema, hydrocephalus, brain herniation, epilepsy, and metastases to other places in several 
areas of the brain, such as the temporal lobe, frontal lobe, parietal lobe [6]. In the case of brain anatomy, MRI 
uses strong magnetic fields, radio waves, and computers to produce more distinct and detailed images of the 
brain and skull structures than other methods. 

Segmentation is one of the methods in digital image processing besides image compression, 
restoration, analysis, and others. The segmentation process of brain tumor MRI images aims to divide the 
area of the MRI image based on the features (color, texture, density, contrast, and other features) [7]. The 
accuracy of the segmentation results can be used for medical analysis processes such as tumor detection 
based on the areas formed. The image segmentation results can also be in the form of an area called a region 
of interest (RoI). The problems faced by MRI image quality for the brain and other organs affect the 
segmentation process to find RoI. Some of them include having bad boundaries. These bad boundaries may 
lead to either inhomogeneity or different intensity ranges, making it challenging to spot the differences 
between pixels in the image [8]. The pixel intensity in the tumor area tends to overlap with pixels in adjacent 
normal tissue. 

The complexity of the MRI image as the structure and morphology of the human brain is very 
complex, brain MRI images tend to have low contrast [9], [10], a relatively large manual segmentation time. 
It requires very high accuracy so that no pixel is missed considering that each pixel contains important 
information, and in some cases, different MRI images show distinct complexity even though they represent 
the same tumor type. In addition, the MRI will produce many images as it is captured from various angles 
(axial, coronal, and sagittal). This can trigger differences of opinion from different experts and the length of 
time it takes to read the entire image. 

Based on the explanations, this study targets to propose a segmenting method for brain MRI images. 
Segmentation of brain MRI images can be served as a crucial early stage to help early detection, increase the 
accuracy and facilitate [11] the process of diagnosing brain tumors. Ultimately, it will help to identify the 
most appropriate type of brain-related treatment. 

The development of the brain MRI image segmentation method begins with low-level techniques, 
such as thresholding [12] and region growing techniques. The thresholding method would group objects by 
comparison with one or more other intensity groups. However, this approach cannot accommodate all the 
pixels of the MRI image. Another method that can be used is the unsupervised learning approach, such as 
clustering [13] and segmentation [14]. The latest segmentation technique used is the deep learning 
technique [15]. Deep learning is the development of artificial neural networks in machine learning. 
Therefore, it can also be called a deep neural network. Deep learning performs the classification of pixels 
during the segmentation process [16]. The variant of the deep learning technique that has been used is the 
convolutional neural network (CNN), which accommodates the formation of a learning layer (network) from 
a set of data automatically. This method classifies each image patch into several classes, such as normal 
tissue, necrosis, edema, and other classes. This classification aims to label the center point of the voxels [17]. 
The hybrid method also developed into wavelet method with machine learning [18], and a combination of 
fully convolutional neural network (FCNN) and conditional random field (CRF) [19]. 

In line with the development of brain MRI image segmentation methods, it is proposed to combine 
two or more methods (hybrid method) to segment images. This study proposes a hybrid method by 
combining the The double density dual-tree complex wavelet transform (DDDTCWT) and CNN with 
additional genetic algorithm (GA) approaches based on its capabilities and mechanisms. 

The combination of these methods will aid each algorithm involved to be more optimal in its 
implementation DDDTCWT algorithm is a technique that combines double density wavelet transform and 
dual-tree wavelet transform [20]. The difference between the two only lies in the number of directions 
(frequency) of the two previous methods so that the results will be more numerous and detailed. The wavelet 
transform involves n levels of decomposition [21] to extract the features contained in the image using the 
filtering function. With more features that can be extracted, it will be easier to carry out the following 
analysis process, including the segmentation process. In this research, DDDTCWT was used to transform the 
input images in detail. Each level of decomposition with DDDTCWT will involve six filtering functions, 
which represent low pass filtering and high pass filtering. This detailed decomposition process will identify 
more feature information contained in the image and increase CNN's recognition and learning process. 
DDDTCWT can also solve the shift-invariance and directional selectivity. Shift invariant occurs when the 
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image is shifted, while the directional selectivity means that more information will be extracted as the 
direction/orientation of the pixel extraction at each level of decomposition. 

CNN algorithm can recognize an image through the learning process (courses or training) with a 
given object. Then, CNN studies all the features in the image through a series of processes, including 
convolution, activation, pooling, and the formation of a fully connected network (FCN) [22]. There were 
additional process options during the formation of FCN, namely the optimization of weights and biases 
involving a GA [23]. Here, we only optimized the first layer of the CNNs with GA and saw its performance. 
GA works as a natural successive that would produce springs, mutation, cross-product, and so on. In many 
cases, the FCN algorithm proved can solve an optimization problem [24]. Then, GA will optimize the 
weights and biases from the first layer of the CNNs layers so that the identification process will be easier and 
faster. In addition, GA will update the parameter values in the training data process by not doing repetitions. 
The principle of GA is to select the gene with the best fitness value so that the weight and bias values 
included in each optimization process (crossovers and mutations) are with optimal values. 


2. LITERATURE REVIEW 

Several studies have been conducted regarding various proposed algorithms for image segmentation. 
These algorithms use diverse approaches, such as thresholding, classification, clustering, and hybrid method. 
Each application of the algorithm gives different accuracy results in the range of 80 - 99%. There are several 
updates in image segmentation methods, from conventional to hybrid methods. Threshold and region-based 
approaches are part of the conventional one. Gaussian mixture model algorithm for segmenting brain MRI 
images is an example of applying a region-based approach [25]. Another study in 2016 that used the 
threshold concept was the magnitude base multilevel method [12]. Moreover, some studies use the concept of 
supervised and unsupervised learning for brain MRI image segmentation using the adaptive K-Means 
Clustering method and the support vector machine (SVM) method [26]. An unsupervised learning method 
that combines fuzzy c-means with k-means clustering has also been studied in India [27] from 2018 to 2020, 
several researchers combined one or more techniques to help identify the presence of brain tumors, such as 
combining a discrete wavelet transform with a probabilistic neural FCN [28] and fully convolutional neural 
network (FCNNs) with CRF [29]. The following Table 1 shows several previous studies with different 
methods. 


Table 1. Segmentation methods used in previous research 


No __ Research Method DSC (%) PPV (%) Sensitivity (%) 
1 Havaei [9] Deep Neural Network 88 89 87 
2 Kamnitsas [29] Fully Connected CRF 89.8 - 89.1 
3 Tustison [29] Concatenad Random Forest 87 85 89 
4 Kwon [29] Combining Generative Model 88 92 84 
5 Pereira [30] Convolutional Neural Network 88 88 89 
6 Zhao [30] FCNNs dan CRF 87 92 83 
Fs Shen [30] Fully Convolution Network 87 85 89 
8 Senthilkumar [31] Back Propagation Neural Network - - 89 
3. METHOD 


3.1. Research framework 

This section contains an explanation of the research framework, the data used, the proposed image 
segmentation method, and the performance evaluation technique of the proposed method. This section also 
describes about input parameter of each algorithm, scenario used for experiment, part of pseudocode for the 
MATLAB code, and CNN architecture proposed. This study has two problem formulations: (1) how to apply 
a combination of three methods (DDDTCWT, CNN, and GA) to segment brain MRI images. It means how to 
combine the workings and input parameters of each method to segment brain MRI images. (2) how is the 
performance of the proposed method to segment the image and identify the grade of tumor or normal tissue. 
The performance was evaluated with four indicators: Dice similarity coefficient (DSC), positive present 
value (PPV), sensitivity, and accuracy as shown in Figure 1. Thus, the purpose of this research is to develop 
the architecture of the proposed hybrid method into a software and evaluate its performance. This study used 
several stages to achieve the goal, including literature reviews, conducting simulations/experiments on the 
developed software, and observing the results of experiments/simulations. 


3.2. Samples 
The images used in this study were grouped into three, namely training images, testing images, and 
ground truth images. The training images were used to study the features contained in all the training images 
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by performing feature extraction. Testing images for examining the image segmentation to produce Rol, and 
the ground truth images were a collection of testing images that have been segmented by the media or other 
methods and were utilized as a comparison with the output image segmented by the proposed method. This 
ground truth image represents the normal image and the tumor image as shown in Figure 2. All MRI images 
were obtained from several hospitals in Indonesia and the Kaggle website. In total, there were 1699 training 
images (1397 MRI images of normal patients and 302 for brain tumor patients) and 913 test images (359 
tumor MRI images and 554 normal MRI images). 913 images from Kaggle sources were used as ground truth 
images to compare the output images of the segmentation with the proposed method. The distribution of 
training images with test images was done randomly by MATLAB with a ratio of 70:30. This study does not 
specifically use the image sharing approach as the amount of data available was adequate. In addition to 
training images and testing images, 


PROBLEM | DATA PROPOSED OBJECTIVE MEASUREMENT 
| METHOD 
¢ Image Training: Building a 
Low in contrast and bad 1,699 of total * DDDICWT hybrid method neers 
boundaries (difficulty images with 1397 ° CNN application . er 1 Score) 
distinguishing normal tissue or normal MRI images * Additional model a: FEN). 
tumor tissue), the quantity of the and 302 brain tumor Genetic @DDTCWT, * Sensitivity 
image, and the potential for MRI images Algorithm CNN, and GA) * Accuracy 
subjectivity in reading the MRI * — Test Image: 913 of to segment brain 
: 7 = total images with MRI images 
ila 359 tumor MRI 
‘ : images and 554 
How to implement a hybrid normal MRI images L . 


method to segment brain MRI 
images 


¢ Ground Truth 
Imagery. A number 
of training images 


Figure 1. Shows the research framework 


Figure 2. Normal and tumor MRI images (Kaggle) 


3.3. Proposed method 

The method proposed in this study is a combination of DDDTCWT, CNN, and GA. DDDTCWT is 
a development of the discrete wavelet transform (DWT) approach with the ability to decompose an image 
into several sub-images up to the nth level. Every decomposition process will produce four sub-bands, which 
are approximation sub-band, diagonal detail sub-band, vertical detail sub-band, and diagonal detail sub-band. 
Image decomposition will be performed on the approximation sub-band because it represents the original 
image. The decomposition results assist in the extraction of detailed information (features) contained in the 
image. The more features obtained, the easier the image segmentation process will be. DDDTCWT improves 
the limitations of DWT in addressing the directional selectivity. DDDTCWT can accommodate a shifting 
image (invariance translation) and more extraction from image pixels (directional selectivity). The input 
parameters used by the DDDTCWT method for the segmentation process are the level of decomposition and 
the wavelet function. In this study, the levels of decomposition used were level 1 to level 5. Meanwhile, the 
wavelet function used three options, including the DDDTCWT filter function (dddtf1), self1 and self2 filter 
functions (for the wavelet transform approach). 
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The same goes for the multi-layer convolution approach in the deep learning method. CNN is 
trained to understand the details of an image better as it works on the principle of multi-layer. Thus, CNN is 
able to learn to identify the different features or details of an image distinctively. The CNN algorithm would 
carry out the learning process of all the features in the training image obtained from the extraction by 
DDDTCWT. After CNN learns the features contained in the training image, 60 FCN scenarios will be 
formed. The input parameters for the CNN algorithm were the convolution layer, activation function, 
pooling, and optimization function. The parameter values of the proposed activation function were rectified 
linear units (RELU) and hyperbolic tangent (TANH) functions, as both are easy to perform and do not cause 
large computations, with the consideration that the data used was adequate. This activation function was 
expected to obtain a convergent state quickly. Table 2 shows a list of the input parameters used from the 
algorithms involved. 


Table 2. Input parameters of each algorithm to form FCN 
No Input Parameter Instance Information 
1 Decomposition Level 1-5 Decomposition level up to 5 according to the size of the 
training image. In addition, a greater level of decomposition 
will take more time. 


2 Wavelet Function DDDTF1 Double Density Dual Tree Filtering 1 (DDDTF1) is the filter 
SELF1 function of the DDDTCWT method, while SELF1 and SELF2 
SELF2 are wavelet functions of other wavelet transforms. 

3 Number of Convolutions According to requirements In each experiment, the number of convolutions can be 


determined as desired. Referring to previous research, the 
number of convolutions can be done more than once. 


4 Convolution Matrix Size 2x2 or 3x3 The proposed matrix size from previous research 
5 Activation Function Relu Function RELU function is the most common used to perform 
Tanh Function activation. 
6 Pooling Technique Max Pooling The second pooling technique is proposed from previous 
Average Pooling research through literature study. 
7 Optimization (GA) Yes (With GA) Yes: To build a neural network, it is necessary to optimize the 
No (NO GA) weight of each layer. 


No: To build a neural network, there is no need to optimize the 
weight of each layer. 


Next, the parameter values for the optimization function (optimization of weights and biases) were 
applied with or without GA, and MAX pooling was used as the default of the CNN algorithm. Max pooling 
is a pooling parameter value that has been widely used in medical image analysis and provides an average 
accuracy of 90-98%. Table 3 shows four FCN scenarios from the training results, which resulting 60 
schemes. 

Before the formation of a neural network, the weight value of each input layer and hidden layer is 
required. During the convolution process, the weight and bias value optimization was performed by applying 
the GA algorithm. GA aids to gain optimal weight and bias values from all iterations (epochs) that occurred. 
Figure 3 is a CNN architecture chart for the training process from the previous method with DDDTCWT. 
Generally, DDDTCWT performs decomposition, filtering, and re-decomposition (reformat of the image). 
The following is the pseudocode for DDDTCWT used in this study. 


Step 1: Decomposition 

Determining the Decomposition Level 

Level 1...4 
Step 2: Filtering 

Wavelet function call 
case 1 
level=1; 
dtcplx=dddtree2 (typetree, grayImage, level, filtername) ; 
dtcplx.cfs{2} = zeros(size(dtcplx.cfs{2})); 
end 
Case 4 
level = 4; 
dtcplx = dddtree2 (typetree, grayImage, level, filtername) ; 
dtcplx.cfs{l} = zeros(size(dtcplx.cfs{1})); 
dtcplx.cfs{2} = zeros(size(dtcplx.cfs{2})); 
dtcplx.cfs{3} = zeros(size(dtcplx.cfs{3})); 
dtcplx.cfs{5} = zeros(size(dtcplx.cfs{5})); 
end 
Step 3: Re-decomposition Image 

dtImage = idddtree2 (dtcplx) ; 
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No SCENARIO 
RELU NO GA 
DDDDTF1 SELF1 SELF2 
Decomposition Level | Decomposition Level | Decomposition Level 1 
Scenario 1 Decomposition Level 2 Decomposition Level 2 Decomposition Level 2 
Decomposition Level 3 Decomposition Level 3 Decomposition Level 3 
Decomposition Level 4 Decomposition Level 4 Decomposition Level 4 
Decomposition Level 5 Decomposition Level 5 Decomposition Level 5 
TAHN NO GA 
DDDDTFI1 SELF1 SELF2 
Decomposition Level | Decomposition Level | Decomposition Level 1 
Scenario 2 Decomposition Level 2 Decomposition Level 2 Decomposition Level 2 
Decomposition Level 3 Decomposition Level 3 Decomposition Level 3 
Decomposition Level 4 Decomposition Level 4 Decomposition Level 4 
Decomposition Level 5 Decomposition Level 5 Decomposition Level 5 
RELU WITH GA 
DDDDTF1 SELF1 SELF2 
Decomposition Level | Decomposition Level | Decomposition Level 1 
Scenario 3 Decomposition Level 2 Decomposition Level 2 Decomposition Level 2 
Decomposition Level 3 Decomposition Level 3 Decomposition Level 3 
Decomposition Level 4 Decomposition Level 4 Decomposition Level 4 
Decomposition Level 5 Decomposition Level 5 Decomposition Level 5 
Scenario 4 TAHN WITH GA 
DDDDTF1 SELF1 SELF2 


Decomposition Level | 
Decomposition Level 2 
Decomposition Level 3 
Decomposition Level 4 
Decomposition Level 5 


Decomposition Level | 
Decomposition Level 2 
Decomposition Level 3 
Decomposition Level 4 
Decomposition Level 5 


Decomposition Level 1 
Decomposition Level 2 
Decomposition Level 3 
Decomposition Level 4 
Decomposition Level 5 
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The performance of the proposed hybrid method was evaluated with four indicators, namely 
DSC [32] or F-1 score, PPV, sensitivity, and accuracy. The first three values measure the results of the 
comparison between the segmented image by the proposed method and the ground truth image, while the last 
indicator for the ability of the algorithm to detect the class (tumor/normal) of the segmented image. These 
four values were obtained through the values of true positive (TP), false positive (FP), true negative (TN), 
and false negative (FN) in the confusion matrix [33]. The formula to measure the four indicator values are, 


2xTP 
DSC = (TP+FN)+(TP+FP) 
Pp = 
(TP+FP) 
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Sensitivity = ai 
y= (TP+FN) 
TP+TN 
Accuracy = ——¥_——_ 
(TP+FP+FN+TN) 
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(3) 


(4) 


Where: 

DSC = Dice Similarity Coefficient 
PPV = Present Positive Value 
Sensitivity = Sensitivity 

Accuracy = Accuracy 

TP = True Positive 

TN = True Negative 

FP = False Positive 

FN = False Negative 


4. RESULTS AND DISCUSSION 

The output of this research is a software of the proposed hybrid method architecture and its 
performance evaluation. MATLAB apps were used to segment brain MRI images as shown in Figure 3 with 
the proposed hybrid method as it provides the toolbox needed and is easy to use, such as wavelet transform, 
machine learning, and deep learning. Furthermore, MATLAB works on a matrix-based, making it easier to 
learn how the compiled program works. Previous studies had published output in an architectural design for 
the combined application of these three methods [34]. The hybrid method architecture was developed 
through three major processes, namely training, testing, and evaluation. The training process aims to conduct 
learning on all training images based on the contained features. Concluding on how each method works, the 
input parameters involved, and the value of each input parameter, this study proposes four scenarios that 
produce 60 FCN combinations. The input parameters for the DDDTCWT algorithm are the level of 
decomposition and the filtering function used. The proposed decomposition levels ranged from | to 5 with 
three filtering functions, namely the filter function DDDTF1, SELF1, and SELF2. DDDTF1 is a wavelet 
filter function recommended for DDDTCWT, while SELF1 and SELF2 are the default filter functions of the 
wavelet transform and are available in the MATLAB library. FCN is a representation of the learning 
outcomes of the proposed method on the features contained in all training images. 

These 60 FCN combinations were then tested on 913 images, yielding the performance value of the 
proposed hybrid method. Tables 4 to 7 respectively show the performance of the methods according to the 
combination in Table 1. Aided by the formed FCN, each test image will be segmented one by one while still 
entering the input parameters according to the FCN used. In other words, this process goals to test the 
learning outcomes (FCN) in recognizing the features and can segment or form Rol from each of the testing 
images. Figure 4 is the graphical user interface (GUI) of the proposed method in MATLAB. Based on 
Figure 4, the left side is the part for entering the input parameter values from the DDDTCWT algorithm and 
calling the FCN that has been generated from the training process. It also presents the performance of the 
proposed method in four indicators, namely DSC, PPV, sensitivity, and accuracy. The right side respectively 
represents the test image used, the segmented image using the DDDTCWT algorithm, Rol from the test 
image of the hybrid method, RoI from the ground truth image, and the segmented image by the hybrid 
method. 


Table 4. Performance of scenario | 


No. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
RELU Not GA 
Measurement dddtf1 selfl self2 
LDI1 LD2 LD3 LD4 LDS LDI LD2 LD3 LD4 \LDS5 LDI LD2- LD3 LD4~ LD5 
98.1 98.5 98.6 98.7 98.7 97.8 98.2 98.5 98.7 98.5 98.2 98.5 98.6 98.3 98.4 
DSC (%) 9 3 1 7 5 7 9 2 0 4 9 ei 3 2 3 
97.5 98.0 98.0 98.3 98.2 96.8 97.6 97.9 98.3 97.9 97.7 98.1 98.2 97.7 97.5 
PPV (%) i 2 2 7 Z 2, 6 4 9 3 0 8 e 1 5 
Sensitivity 98.9 99.1 99.2 99.2 99.3 99.0 99.0 99.1 99.0 99.2 98.9 990 99.1 990 99.3 
(%) 6 1 8 4 Z 2, 1 6 7 2 he 2 0 7 9 
96.5 97.1 97.3 97.6 97.6 95.9 96.7 97.1 97.5 97.1 96.7 97.2 97.3 96.7 97.0 
Accuracy (%) 4 7 4 5 0 3 2 af 1 9 1 4 if 9 0 
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Table 5. Performance of scenario 2 


No. 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 


TANH Not GA 
Measurement dddtf1 self self2 
LDI1 LD2 LD3 LD4 LDS LDI LD2 LD3 LD4 LDS LDI LD2- LD3 LD4~ LD5 
97.1 98.2 98.0 97.9 98.4 98.2 97.8 98.1 98.4 980 98.1 98.0 98.0 97.8 97.8 


DSC (%) 8 1 4 9 3 9 3 7 3 8 a) 3 7 2 1 
95.4 97.2 96.8 96.8 97.5 97.6 96.6 97.1 97.8 96.9 O73 97.0 97.0 96.6 96.4 
PPV (%) 9 8 5 3 6 7 1 7 7 0 7 1 1 2 6 
Sensitivity 99.0 99.2 99.3 99.2 99.3 98.9 99.2 99.2 99.0 99.3 99.0 99.1 99.2 99.1 99.3 
(%) 9 4 7 9 9 9 0 6 8 6 3 8 3 6 2 
94.6 96.5 96.2 96.1 97.0 96.7 95.8 96.4 97.0 96.3 96.4 96.2 96.3 95.8 95.8 
Accuracy (%) 9 8 8 7 0 1 7 9 1 4 8 5 2 7 5 


Table 6. Performance of scenario 3 


No. 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 
RELU WITH GA 
Measurement _ dddtf1 selfl self2 
LD1 LD2 LD3 LD4 LDS LDI LD2- LD3~ LD4 LDS LDI LD2 LD3 = LD4~ LD5 
98.3 98.8 98.4 986 98.5 98.2 98.2 98.6 98.0 984 98.5 98.7 98.3 98.6 
DSC (%) 5 4 5 3 Wi 8 8 4 99.42 2 6 6 1 9 9 
97.8 986 976 982 979 976 97.5 98.2 100.0 97.0 98.0 98.1 98.3 97.8 98.3 
PPV (%) 6 2 6 0 5 9 6 7 0 5 3 9 9 3 9 
Sensitivity 98.9 99.1 99.3 99.1 99.2 989 990 99.0 990 989 990 990 99.0 99.0 
(%) 4 1 2 3 5 6 8 8 98.86 9 5 0 7 2 6 
Accuracy 96.8 97.7 97.0 97.33 97.2 96.7 96.7 97.3 96.1 97.0 97.2 97.5 96.9 974 
(%) 5 7 3 7 5 0 0 8 98.86 9 4 3 1 1 9 


Table 7. Performance of scenario 4 


No. 46 47 48 49 50 51 52 53 54 55 56 af 58 59 60 
TANH WITH GA 
Measurement __dddtf1 selfl self2 
LDI LD2 LD3_ LD4 LD5 LDI  LD2 LD3) LD4 LDS5 LDI LD2 LD3 LD4_~ LDS5 
975 984 98.4 98.6 97.5 98.1 98.3 98.1 96.5 98.2 986 980 95.9 
DSC (%) 5 4 i 99.42 1 a 99.42 2 2 2 I ee 3 2 3 
96.2 97.7 97.5 100.0 97.9 96.2 100.0 97.0 97.5 97.3 94.1 974 98.1 96.9 93.0 
PPV (%) 5 3 7 0 5 4 0 9 3 2 a 1 4 8 4 
Sensitivity 990 99.2 99.3 99.3 99.0 99.3 99.1 99.0 99.2 99.1 99.1 99.1 99.1 
(%) 1 5 4 98.86 4 1 98.86 0 8 2 3 w2 6 8 6 
Accuracy 95.3 970 96.9 97.3 95.3 964 967 964 935 965 97.3 96.2 92.3 
(%) 5 1 7 98.86 4 8 98.86 5 8 0 3 9 7 3 2 


The next output is from measuring the performance of the proposed method for segmenting and 
detecting the testing images of tumor or normal. The measurement indicators used were DSC, PPV, 
sensitivity, and accuracy. Based on the first experimental scenario with the FCN formation scheme of RELU 
activation function without involving GA, the highest DSC and accuracy values were in DDDTF1 wavelet 
function and decomposition level 4 with 98.7% and 97.65%, respectively. However, the highest sensitivity 
value (98.39%) was in SELF1 with decomposition level 5. The utmost PPV value was obtained when using 
the SELF] filter function with decomposition level 4 with 98.39% as shown in Table 4. Then, the 
performance of the second scenario, the TANH activation function without GA, shows the highest DSC, 
PPV, and accuracy values with the fourth decomposition level and wavelet function SELF1, while the 
highest sensitivity value was in the DDDTF1 of the fifth decomposition level as shown in Table 5. Next, the 
third scenario involves the RELU activation and the GA optimization function in the FCN formation. Table 6 
shows that SELF! with the fourth decomposition level had the highest DSC, PPV, and accuracy values, 
while the DDDTF1 with the third decomposition level had the highest sensitivity value. The last scenario 
involves TANH activation function and the GA optimization had two schemes with the highest DSC, PPV, 
and accuracy values, namely DDDTF1 with decomposition level 4 and SELF1 with decomposition level 2. 
Sedangkan nilai akurasi tertinggi ketika menggunakan fungsi filter DDDTF1 dan level dekomposisi 3 as 
shown in Table 7. 

Overall, only scenarios 3 and 4 have PPV values reaching 100%. It proves that the integration of the 
GA optimization function affects the PPV value, and the scheme without GA shows that the largest PPV 
value is always at the fourth decomposition level. The highest DSC value also occurs in those scenarios when 
applying both RELU and TAHN activation functions with GA optimization. Likewise, the highest sensitivity 
value falls in the scheme of THAN activation function with and without GA. In detail, the highest sensitivity 
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value lies in the DDDTF1 with the third decomposition level. Lastly, the best accuracy values of 98.86% lie 
in scenarios three and four, precisely in the SELFl with decomposition level 2 (scenario 3), DDDTF1 with 
decomposition level 4 (scenario 3), and DDDDTF1 filter function with decomposition level 4 (scenario 4). 

Based on its performance, all measurement indicators of the proposed method had scores ranging 
from 95% to 100%. This shows that the proposed hybrid method provides a better contribution than previous 
studies to perform image segmentation. This also proves the ability of each method used, CDDDTWT as a 
variant of the wavelet transform approach, to extract all the features contained in all training images. 
Moreover, CNN can learn distinctively all features in the images and can produce optimal weight and bias 
values when integrated with GA. 
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Figure 4. Application of the proposed testing method with the MATLAB tool 


5. CONCLUSION 

The combination of DDDTCWT, CNN, and GA can be applied to segment brain MRI images. The 
combination of the mechanisms and capabilities of the three methods shows the success of segmenting brain 
MRI images with excellent assessment indicators. All combinations of the proposed hybrid system in this 
study showed more than 95% in all parameters measured (DSC, PPV, sensitivity, and accuracy) using 913 
test images. The top 3 combinations were (1) RELU with GA, filter type of DDDTCWT in SELF1 with 
decomposition level 4, (2) TANH with GA, filter type of DDDTCWT in DDDTF1 with decomposition level 
4, and (3) TANH with GA, filter type of DDDTCWT in SELF1 with level decomposition 2. The smallest 
value found at TANH with GA, filter type of DDDTCWT in SELF2 with level decomposition 1. Adding the 
GA in the system improved the system measurements on average, although the time consumed by the system 
to apply GA is doubled. Further improvements of this hybrid method with another technic and CNN training 
are required for a more applicable technique. In addition, this study also suggests increasing the number of 
data sets, both training and testing data. Furthermore, if the number of data sets obtained is quite large, it 
needs the cross-validation technique in the distribution of the image, given that cross-validation can increase 
the model's performance. 
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